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This paper describes the fundamental principles of launch vehicle flight control analysis and design. In 
particular, the classical concept of “drift-minimum” and “load-minimum” control principles is re-examined 
and its performance and stability robustness with respect to modeling uncertainties and a gimbal angle con- 
straint is discussed. It is shown that an additional feedback of angle-of-attack or lateral acceleration can 
significantly improve the overall performance and robustness, especially in the presence of unexpected large 
wind disturbance. Non-minimum-phase structural filtering of “unstably interacting” bending modes of large 
flexible launch vehicles is also shown to be effective and robust. 


I. Introduction 

Note to Session Organizer/Reviewers: This draft manuscript summarizes very preliminary results obtained dur- 
ing an early phase of a project for the launch vehicle flight control systems analysis and design as applied to Ares-I 
Crew Launch Vehicle. During the next several months, a more detailed, rigorous study will be conducted in the areas 
of drift-minimum vs load-minimum control, flexible-body stabilization and analysis, gain scheduling vs. adaptive 
control, etc. A companion paper on dynamic modeling of large flexible launch vehicles is also being submitted to this 
Space Exploration and Transportation GNC session. 


II. Rigid-Body Control Analysis 

Consider a simplified linear dynamical model of a launch vehicle [15], as illustrated in Fig. 2, as follows: 

0 = M a a + Mgd (1) 

Z = - —i 9-—a+—S (2) 

m m m 

z 

& = 0 + — + a w (3) 

F = T a +T~ D (4) 

where 9 is the pitch attitude, a the angle of attack, Z the inertial Z-axis drift position of the center-of-mass, Z the 
inertial drift velocity, m the vehicle mass, T a the ungimbaled sustainer thrust, T the gimbaled thrust, N = N a a the 
aerodynamic normal (lift) force acting on the center-of-pressure, D the aerodynamic axial (drag) force, F the total 
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Figure 1. Comparison of Space Shuttle, Ares I, Ares V, and Saturn V Launch Vehicles [1J. 


x-axis force, S the gimbal deflection angle, V the vehicle velocity, a w = V w jV the wind-induced angle of attack, V w 
the wind disturbance velocity, and 


M a =x cp N a /Iy (5) 

M S = X C gT/Iy (6) 

N a = i pV 2 SC Na (7) 

where I y is the pitch moment of inertia. For an effective thrust vector control of a launch vehicle, we need 

M s S 

max >M a 

^ma x (8) 


where <J max is the gimbal angle constraint and a max is the maximum wind-induced angle of attack. 
The open-loop transfer functions from the control input 8(s) can then be obtained as 
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(12) 


Consequently, the 4th-order system described by Eq. (1) - (3) is completely controllable by <5 and is observable by Z\ 
however, the system is not observable by 6 and a. 

In 1959, Hoelkner introduced the “drift- minimum” and “load-minimum” control concepts as applied to the launch 
vehicle flight control system [6]. The concepts have been further investigated in [7-14]. Basically, Hoelkner’s con- 
troller utilizes a full-state feedback control of the form 


8 = -K±e - K 2 9 - K 3 a 


2 of it 


(13) 
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Figure 2. Simplified pitch-axis model of a launch vehicle. 


for a 3rd-order dynamical model of the form 
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This 3rd-order system is observable by 9 or a. The feedback gains are to be properly selected to minimize the lateral 
drift velocity Z = V(a — 9 — a w ) or the bending moment caused by the angle of attack. Note that 


— =7 = a-9-a w (15) 

where 7 is often called the flight-path angle. 

Instead of measuring the angle-of-attack, we may employ a body-mounted accelerometer, as illustrated in Fig. 2, 
as follows: 


S = - K x 9 - I< 2 9 + K a z m 

= - J\ | 9 - Ko9 + K„ (-—a +-6+ —fi) 

\ m in m ) 

r ■■ N T 

= - K[ 9 - K 2 9 + K a — 9 - K a — -a + I< a -S 

m m in 

Because the resulting effect of z m feedback is basically the same as the a feedback, we consider here only the control 
logic described by Eq. (13). 

Substituting Eq. (13) into Eq. (1) - (2) or Eq. (14), we obtain the closed-loop transfer function from the wind 
disturbance a w (s) to the drift velocity Z(s ) as 

Z _ A 2 s 2 + A\s + A a 

a^V ~~ s 3 + B 2 s 2 + B lS + B 0 { ’ 
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For a unit-step wind disturbance of a w (s) = 1/ s, the steady-state value of Z can be found as 

Z ss _ — (A 2 S 2 + A\s + A 0 ) — —A a _ 1 

V s^o s 3 + B 2 S 2 + Bis + B 0 B 0 


(17) 


The launch vehicle drifts along the wind direction with Z ss = — V v , and also with 9 = 9 = a = 5 = 0ast — + 00 . 
It is interesting to notice that the steady-state drift velocity (or the flight path angle) is independent of feedback gains 
provided an asymptotically stable closed-loop system with B 0 0. 

If we choose the control gains such that B 0 = 0 (i.e., one of the closed-loop system roots is placed at s = 0), the 
steady-state value of Z becomes 


Z S S _ y —{A 2 S+A 1 ) _ ~A\ _ —1 

~V - '™o S 2 + B 2 s + B 1 ~ ~B/~ ~~ 1 + c 


(18) 


where 

mVlMsilu +K 3 )~ M a ] 
M a K 2 T+M s N a /T 


For a stable closed-loop system with Mg(K\ + K 3 ) — M a > 0, we have C > 1 and 


(19) 


\Zss\ < v w 


( 20 ) 


when B 0 = 0. The drift-minimum condition, B a = 0, can be rewritten as 


MgK 3 — M a N a f x cp \ 

M S K ! F V x cg ) 

Consider the following closed-loop transfer functions: 

(x s(s 2 -(- MgK 2 s MfiNi'j 
cx w s^ 4- B 2 s 2 + B\S 4- B 0 


( 21 ) 


( 22 ) 


= s(K 3 s 2 + M a K 2 s + Mplu) 
a w s 3 + B 2 s 2 + Bis + B 0 

For a unit-step wind disturbance of a w (s ) = 1/s, we have a = 6 = 0 as t— > 00 . However, for a unit-ramp wind 
disturbance of a w (s ) = 1/s 2 , we have 


lim a(t ) = M$Ki 

t — >00 

lim S(t) = M a K 1 

t—> OO 

Consequently, the bending moment induced by a and S can be minimized by choosing K\ = 0, which is the “load- 
minimum” condition introduced by Hoelkner [6], The closed-loop system with K 1 = 0 is unstable because 
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mV 


(. M S K 3 - M a ) < 0 


(24) 
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However, the load-minimum control for short durations has been known to be acceptable provided a deviation from 
the nominal flight trajectory is permissible. 

A set of full-state feedback control gains, ( K \ , A 2 , A 3 ) , can be found by using a pole-placement approach or the 
linear-quadratic-regulator (LQR) control method [21-22], as follows: 

/»oo 

min / [x t Qx + S 2 )dt (25) 

s Jo 

subject to x = Ax + BS and <5 = —Kx where x = [9 9 a] T and A = [K 1 K- 2 A 3 ]. 

III. Rigid-Body Control Example 

Consider a launch vehicle control design example discussed by Greensite in [15]. Its basic parameters are given as 
in [15] 


Iy = 2.43A6 slug-ft 2 , 
F = 375, 000 lb, 

V = 1320 ft/sec, 

N a = 240,000 Ib/rad, 


to = 5830 slug, 
x cp = 38 ft, 

V w = 132 ft/sec, 
M a = 3.75 s- 2 , 


T = 341,000 lb 
x cg = 32.3 ft 
a w = 5.73 deg 
M s = 4.54 s " 2 


The open-loop poles of this example vehicle are: -1.9767, 0.0488, 1.8967 


(26) 


Note that the wind-induced angle of attack of 5.73 deg considered for this example in [15] is somewhat unrealistic 
because it will require a maximum gimbal deflection angle of 

c Ma . 7 , , 

Umax ^ , r Oi w — 4.73 deg 
Ms 

Most practical thrust vector control systems have a maximum gimbal angle constraint of about ±5 deg. In this paper, 
we also assume a second-order gimbal actuator dynamics of the form 


5(s) = 


U)Z 


S 2 ' h 2 0J n s -\- i Uj 


;5 c (a) 


(27) 


where £ = 1 and uj n = 50 rad/s. 


Table 1. Summary of rigid-body control analysis and design 


Case No. 

Controller Type 

Feedback Gains ( A’i . K 2 . A 3 ) 

Closed-Loop Poles 

1 

(9, 0)-Feedback Control [15] 

( 2 , 0 . 8 , 0 ) 

-1.7488±1.3934j, -0.1596 

2 

Drift-Minimum Control [15] 

(2,0.8,3.614) 

-1.9087±4.2774j, 0.0 

3 

Load-Minimum Control [15] 

(0,0.8, 3.614) 

-1.9323±3.0533j, 0.0471 

4 

LQR Control (Q = 0) 

(0.6852,0.8491,0.9542) 

-1.9767,-1.8967,-0.0488 

5 

Drift-Minimum Control 

(0.3220,0.8352, 1.2765) 

-1.9767,-1.8967,0.0 

6 

Load-Minimum Control 

(0, 0.8352, 1.2765) 

-3.1323,-0.7816, 0.0405 


IV. Flexible-Body Control Analysis 

More detailed control and stability analysis results for Figs. 10 and 1 1 will be included in the final manuscript. 
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Figure 3. (0, 0) -feedback control (Case 1). 
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Figure 4. Drift-minimum control (Case 2). 


















Mode Disp Slope 


O) 

133 „ 
"D 4 


2 

I 

~ °P 

D) 

03 C 

T3 -5 ' 


- 



1 

1 

- 

0 

1 

2 

3 

4 

5 


-10 L 


05 

03 n 
T3 0 



1 

0.5 
0 

-0.5 

-1 

500 1000 1500 2000 2500 3000 3500 

Sensoi Location Fiom Tip [in] 


-4 

x 10 




t 

Instrument 

Unit 


j^SBSS T*:'T IIIII1. 1 1 


t 

Upper Stage 
Aft Skirt 


t 

First Stage 
Aft Skirt 


Figure 10. Illustrations of dominant bending modes and sensor locations (Ref. 2). 
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Nichols Chart: Flight Time at 60 sec 



Figure 11. Nichols plot for a baseline pitch-axis flight control system (Ref. 2). 
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V. Conclusions 
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